Average magnification effect of clumping of matter 



T.W.B. Kibble 

"i/-^" ■ Blackett Laboratory, Imperial College, London SW7 2AZ, United Kingdom 

o : 

O' kibbleOimperial . ac .uk 

(N ; ^ 

S ! and 
1—5 ; 

. Richard Lieu 

(N : 

^ ' Department of Physics, University of Alabama at Huntsville, Huntsville, AL 35899, U.S.A. 

> ■ 

ly-^ ■ lieur@cspar.uah.edu 

i> : 

(N ■ 

<N ■ ABSTRACT 

■ The aim of this paper is to re-examine the question of the average magnifi- 
Q^' cation in a universe with some inhomogeneously distributed matter. We present 
Q ■ an analytic proof, valid under rather general conditions, including clumps of any 

■ shape and size, and strong lensing, that so long as the clumps are uncorrelated 
• the average reciprocal magnification (in one of several possible senses) is precisely 

■ the same as in a homogeneous universe with equal mean density. From this re- 
^ ■ suit, we also show that a similar statement can be made about one definition 

■ of average direct magnification. We discuss, in the context of observations of 
discrete and extended sources, the physical significance of the various different 
measures of magnification and the circumstances in which they are appropriate. 

Subject headings: gravitational lensing — galaxies, observations — distance scale, 
distances and redshifts — cosmology 



1. Introduction 

There has been considerable debate about the average magnification effect of gravita- 
tional lensing by randomly distributed clumps of matter. Weinberg (13) argued that the 
average magnification produced by randomly distributed masses is exactly the same as that 
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in a homogeneous universe of equal mean (or pre-clumping) density — the magnification 
produced by the clumps is largely cancelled by the Dyer-Roeder effect (4; 5). But his argu- 
ments have been criticized by Ellis et al. (6), who pointed out that they ignore the effects of 
caustics. These authors also introduced an important distinction between two measures of 
distance, which they called 'area distance' and 'angular-size distance', though in fact both 
can be applied to either lengths or areas. Holz & Wald (8) developed a general formalism 
for estimating the probability distribution of magnification, as well as shear and rotation, 
and obtained numerical results for a range of cosmological parameters using Monte Carlo 
simulation of light paths. Claudel (3) studied a number of different examples, and con- 
cluded that to first order small deviations from homogeneity would not change the average 
magnification. On the other hand. Rose (12) gave an analytic argument using a spherically 
symmetric model of the universe with the aim of showing that objects in an inhomogeneous 
universe appear, on average, more magnified than those at the same redshift in a homo- 
geneous universe with the same mean density. This is not in contradiction with Claudel's 
result, because the effect Rose finds is of second order. 

The purpose of this paper is to re-examine this question using a simple and explicit 
analytic approach. We show that under rather general conditions there is at least one 
measure by which the average reciprocal magnification is exactly the same as in an FRW 
universe with the same mean density. When there is strong lensing, the different measures 
of distance diverge. It is easier to deal initially with reciprocal magnification, because it goes 
to zero rather than infinity on the caustics. Later, however, we do consider average direct 
magnification. 

Our starting point is in some respects similar to that of Holz & Wald (8), based on 
using the geodesic deviation equation to follow the paths of light signals back in time. Our 
goal is more restricted, in that we focus only on average magnification, not rotation and 
shear. On the other hand, we are seeking analytic rather than purely numerical estimates, 
so the assumptions we make are slightly more restrictive, though still, we believe, of wide 
applicability. 

Specifically, we assume that in addition to a smooth, homogeneous matter component, 
with density p/^, there is another component comprising widely separated, slow moving and 
randomly distributed mass clumps (say galaxies, groups, or clusters). For simplicity, we 
suppose initially that each clump has the same mass M. But it is easy to generalize the 
discussion to include a distribution of masses, even an evolving one. 

Holz and Wald assumed that the universe can be described by a 'Newtonianly perturbed 
Friedmann-Robcrtson- Walker (FRW) universe' (7), i.e. the metric is an FRW metric with 
the time and space parts multiplied respectively by (1-1-20) and (1—20), where the convention 
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of c = 1 is adopted here and henceforth. With various assumptions on and the matter 
distribution, they showed that (j) obeys a Poisson equation, with 8p = p — p on the right hand 
side, where p is the density and p the density of the corresponding FRW universe defined by 
setting = 0. They argue that to determine the way a hght signal propagates it is sufficient 
to look explicitly only at the gravitational potential of nearby clumps. 

According to our assumptions, the density perturbation would comprise two contri- 
butions, a spatially uniform negative background ph — p, and an occasional large positive 
contribution from one of the clumps. For most of its journey, a light signal will be travelling 
through a uniform background, but when it docs pass near a clump the effects will be much 
larger. Under these conditions it is reasonable to assume that we can deal with the effects of 
the clumps individually. We assume that the clumps are small and slow-moving enough that 
the gravitational effect of each one may be treated in a Newtonian approximation, with a 
time-independent Newtonian potential Moreover, we use the 'plane lens approximation', 
that is, we compute the angle of deviation due to the clump by integrating the gradient 
of the potential along the undcviated light path, and assume that the deviation effectively 
occurs at the central plane. As pointed out by Metcalf & Silk (11), this induces a small 
error because the true light path passes closer to the center. However the discrepancy in the 
minimum distance from the center is very small, of order the Schwarzschild radius of the 
clump. Hence the error is tiny and consistently negligible in the Newtonian approximation. 
Finally, we also assume that far from the clumps there is no appreciable source of shear, so 
that the Weyl tensor vanishes. Of course, no such assumption is made about the field near 
each clump. 

One criticism that might be made is of our assumption that the clumps are well separated 
and randomly distributed. This does not mean however that only one clump can significantly 
affect a light signal at any time (though that may often be true) , but rather that the effects of 
different clumps are purely additive. This seems to us generally a good approximation. The 
most serious objection would probably be to the assumption that the clumps are uncorrelated. 
Such correlations may invalidate the assumption that there is no source of shear far from 
the clumps. Even in such cases, the effect on average magnification should be small, since 
according to the Raychaudhuri equation the effect of shear on expansion is of second order. 
These correlations might also be thought to call in question the vahdity of the plane-lens 
approximation, but this would be true only if the clumps are correlated in such a way that 
the deviated light paths sample a significantly different environment. Given the extremely 
small error in the deviation angle (typically of order ■0^) , this seems very unlikely. 

It is important to note that the 'average magnification' for a given red-shift can mean 
several different things. In the strong-lensing case, when caustics are present, imaged areas 
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fold back on themselves. In one sense, the magnification is negative in the region beyond 
the caustic, because images are reversed. In the distinction made by Ellis et al. (6), in 
computing the angular-size distance, these regions are indeed counted negatively, whereas 
the area distance is concerned with the total area, including all the folds; in that case, every 
contribution is taken positively. 

There is another important distinction to be made. We may choose at random one of 
the sources at red-shift z, or we may choose a random direction in the sky and look for 
sources there. These are not the same; the choices are differently weighted. If one part of 
the sky is more magnified, or at closer angular-size distance, the corresponding area of the 
constant-^; surface will be smaller, so fewer sources are likely to be found there. In other 
words, choosing a source at random will give on average a smaller magnification or larger 
angular-size distance. 

Which of these definitions is appropriate depends on what we choose to look at, and 
what questions we want to ask. We shall return to the question of which definition to use in 
various circumstances in a later section. 

Let us concentrate for the moment on the random-direction averaging. The question we 
wish to address is this: How is the average magnification affected by whether the matter is 
clumped rather than smoothly distributed? We do this by examining the geodesic deviation 
equation in the presence of clumps. 

One other preliminary point should be made. What we arc interested in obscrvationally 
is the average magnification of sources at a given red-shift z. But what we actually calculate 
is the average of sources at the same affine distance A (along the backward null geodesies 
from the present), which is not exactly the same thing. We argue however that the difference 
is undetectably small. The effect of passing near a clump of mass M affects the relationship 
between z and A in much the same way as the conventional gravitational time delay. Thus 
the difference in z for fixed A is of order HqGM times a logarithmic factor, which is negligible 
under any reasonable conditions. 



The Robertson-Walker line element for an open universe, with k — — \k\ and c = 1, is 



2. Null geodesies 




(1) 
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or equivalently, with t — J dt/a{t) and r = |A;| smh(|A;|^/^x), 

ds^ = (r) I dr^ - dx' - ^ sinh^ {\k\ {d-d^ + sin^ dc^^-, | _ 
Of course in the flat-space hmit, \k\ ^ 0, r and x become identical. 



The Priedmann equation is 
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where p„i is the density of matter (assumed pressureless). Consequently, the relation between 
the Hubble parameter H and the red-shift, z = ao/a{t) — 1, is H = HqE{z), where 

E''{z) = Qmil + Zf +{l-nm- ^a){1 + Zf + ^A, (4) 

in which as usual VL^ = STiGpmo/^HQ and fl\ = A/3Hq. 

We shall consider backward null geodesies from the origin at the present time Iq, with 
affine parameter A normalized so that t{0) — —1, where the dot denotes a derivative with 
respect to A. Then 

ao 

We now assume that in addition to a uniform distribution of matter there are random 
clumps present. Specifically, the matter density parameter may be written = flh+flg, 
where flh represents a homogeneous distribution and fig a random distribution of widely 
separated clumps, each of mass M. (It is easy to generalize the discussion to a distribution 
of masses, or even to allow for a distribution changing with cosmic time.) 

Consider a fiducial backward null geodesic, and a second neighboring null geodesic from 
the same point. We choose a vierbein e(^) at the origin, with e(o) in the t direction, and 
6(3) — 6(0) tangent to the fiducial geodesic. Then we parallel-propagate the vierbein along 
this geodesic, and introduce transverse coordinates I — (/^,i^), such that the transverse 
separation between the geodesies at affine distance A is 

Sx'^ = rel^{X). (6) 

(Summation over a = 1, 2 is implied.) The rate of change of this separation is governed by 
the geodesic deviation equation, 

r = R"^f,.uVl^, (7) 
where — is the null tangent vector to the fiducial geodesic. 
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It is convenient to define the 2x2 matrix ^(A) = ^^"^(A)^ by 

r(A) = ^(A).r(o), (8) 

where /(O) defines the initial (small) angular deviation from the fiducial geodesic. Then A 
also satisfies the geodesic deviation equation (7). It is also useful to introduce the quantity 

A(A) = det^(A), (9) 

which is proportional to the reciprocal magnification of a small source. In the weak-lensing 
case, A is always positive. Then the angular-size (or area) distance D is given hy D = \I~A.. 
However, caustics are lines of infinite magnification, where yl = 0. In the regions beyond, 
where A < 0, images are reversed. More generally. 



D{\) = Vl^l- (10) 

For a single image of a point source, there is no need to distinguish angular-size and area 
distances. The distinction for large sources is discussed in Section 5. 

So long as we are well away from the clumps of matter, it is reasonable to assume that 
the Weyl tensor is very small, z.e., there is no source of shear. We make no such assumption 
about the effect of individual clumps. (Each clump of course introduces a source of shear in 
its vicinity.) Then the coefficient on the right hand side of the geodesic deviation equation 
involves only the Ricci tensor. Using Einstein's equations, this may be expressed in terms 
of the density ph of the smooth component, or equivalently Q/j. Thus, along light paths that 
come near to none of the clumps, we find for A the equation 

A=-\El^n{^ + zfA. (11) 

Here a factor (1 + z)^ comes from the 2;-dependence of p/i, while (1-1-2;)^ derives from the 
u^v!" factor. Clearly, away from the clumps, A is the solution of this equation (11) with 
initial conditions 

^(0) = 0, ^(0) = 1. (12) 

Thus for these paths, we find ^(A) = D/i(A)l, where Dh{X) is the Dyer-Roeder distance, 
the solution of the equation 

Dn^-^H',Qn{l + zfDn, (13) 

with initial conditions 

Dh{0) = 0, Dh{0) = 1. (14) 
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It will also be useful to define the corresponding angular-diameter distance -D(A) for a 
homogeneous FRW universe with density p, which is in fact given by 

D{X) = a{X)r{X). (15) 

If the relative magnification of a point source, compared to that of a source at the 
same red-shift in the FRW universe is defined to be negative when the image is reversed, 
then 

The magnification in the more usual sense, which is always positive, is then 

i„i=iw=m (17) 



3. Effect of a clump 

We now consider the effect of a clump near our null geodesies at affine distance A. We 
suppose that the center of the clump is at transverse position —x from the fiducial geodesic. 
The effect of the clump will be to bend the geodesic through an angle ip. With the center of 

— * 

the clump as origin, let us introduce Euclidean coordinates {x, Q . Then -0 is given in terms 
of the Newtonian potential $(x, () of the clump by 

/oo 
dCV$(a-,C), (18) 
-oo 

where V is the two-dimensional gradient operator. Note that if p{x, Q is the density distri- 
bution in the clump, then, by virtue of the Poisson equation, 

V-'0(f) =47rG'(7(f), (19) 

with 

/oo 
dCp{x,C). (20) 
■oo 

The angular deviation of a neighboring geodesic at transverse displacement /, relative 
to that of the fiducial geodesic, will be 



■0(x + I) — i^ix) — lC{x) • z. 



(21) 
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where the 2x2 matrix K, is given by 

]C{x) = Vif{x). (22) 

Note that by the form of (18) JC is symmetric. 

Now from (5) and (8), the angular deviation of neighboring geodesies at position A is 

i = ^ = M2). (23) 

ax 1 + z 1 + z ^ ^ 

Hence A changes sharply as the geodesic passes the clump by the amount 

6A^ {l + z)}C{x)-A{X). (24) 

It is also possible to find a similar expression for 5 A. Since 

A^Atr{A-^A), (25) 

it follows that 

SA^ A tr{A-^S A) ^{1 + z) A tr}C{x). (26) 

Our aim now is to compute the expected value of the change in ^ in a small interval 
dX due to the effect of a random distribution of clumps. 

Since the clumps are of mass M and constitute a fraction Qg of the critical density, their 
number density (assuming no evolution) is clearly 

'^-Ml + ^f, no ^5^. (27) 

The probability of finding a clump with center at the position (A, —x) within small ranges 
dX, (Px is 

nQ{l + zYdXd^x, (28) 

where the extra factor of z) occurs because adx = (1 + z)dX. 

To find the equation we are looking for, let us ask how {A) changes under a small change 
dX in A. Of course, due to the smooth background, according to (11), there is a change 

[d(^(A))], = —Hl^nil + ;.)^(^(A))ciA. (29) 
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Now what is the average effect of the clumps? The probabihty that a clump is present 
is given by (28). Since we have assumed the clumps are uncorrelated, this probability is of 
course independent of the previous history, so by (24) the change due to clumps is 

[d{A{X))]g = no(l + z)^dX J d^xlC{x) ■ {AiX}). (30) 

The next step is to compute the x integral. Using (22), we have 

I d'^xlC{x) = j d^xV^ix) = ^dlntfix), (31) 

where n — x/\x\, and the integral is around a large circle. Now for large x we have asymp- 
totically 

if{x) ~ -4GM^, (32) 

whence 

J d^xlC{x) = -AttGMI. (33) 

For later use, we need one other result, concerning the integral of the determinant of /C. 
If ^ is a unit vector in the direction of the fiducial null ray, then clearly 

C det /C = VV'i X = V X (V'l VV'2) . (34) 

Hence 

d'^x det/C = i dl ■ (35) 



where the integral is again around a large circle. Clearly this vanishes in the limit, whence 

d'^x det /C = 0. (36) 



Now it follows from (30), (27) and (33) that 

[diAm = -^^^0^^.(1 + z)'{A{X))dX. (37) 

Remarkably enough, this is of precisely the same form as (29), so, combining the two, 
we find 

{A{X))^~H',nr^{l + zr{A{X)). (38) 

In other words, (.4.(A)) is exactly the same as it would be in a homogeneous universe of equal 
mean density, namely 

{A{X)) ^ D{X)1, (39) 
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where D is given by (15). 

It is not possible to use (26) directly to find a similar equation for the mean value of the 
quantity A defined by (9), because even in the homogeneous background A does not satisfy 
a simple homogeneous differential equation hke (11). Nevertheless, we can also show that 



The general proof of this result is rather long and is therefore relegated to an appendix. 
However, in the next section we explain why it is true in the special case where a typical 
light path encounters no more than one clump. 

This result, (40), shows that in a clumpy universe the average reciprocal magnification 
(1///), where /i is defined by (16), is exactly the same as in a homogeneous universe of equal 
mean density. 

It is remarkable that, despite the quadratic relationship between A and A, the averages 
of both, (A) and (A), can be expressed in terms of this one function D. It is very important 
to note, however, that D is not the average (D) of the angular-size or area distance of point 
sources — and nor is D"^ equal to (D^) — because according to (10), (D^) = (|^|)- As we 
shall see, for a large source, D is approximately the 'angular-size distance' of Ellis et al. (6), 
not the 'area distance'. 



To clarify the reason for the important result (40) we shall in this section speciahze to 
the single-clump case, where the effect of multiple encounters is assumed to be neghgible. 
This will make it easier to estimate the size of the differences. 

Let us begin by finding a more explicit expression for (^(A^)) for sources at affine dis- 
tance As in the case where the probability of encountering more than one clump is negligible. 
It will be useful to extend the definition of the angular-diameter distance in the 'empty' 
regions, given by (13) and (14). We define Dfi{\i,X) to be the angular-diameter distance 
between Ai and A in the smooth background between the clumps, i.e., the solution of (13) 
with initial conditions 



(^(A)) = [D{X)r. 



(40) 



4. The single-clump case 



d_ 
dX 



Dh{X,,X) 



A=Ai 



(41) 



(The factor of (1 -|- Zi) appears because ax = 1 -I- z.) 
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If there are no clumps near the hght path, then of course ^(A^) = D^(As)l. If there is a 
single clump at A and transverse position —x, then its deviating effect is added, so one finds 



or equivalently 
where 



^(A,) = DhiK)! + DhiX, Xs)}Cix)DhiX), 
Dh{X)Dh{X,Xs) 



DhiXs) 



It will be useful to define the quantity 

J = det[l + L(A,A,)/C(f)], 

so that 

A{X,) = DliXs)J, 



(42) 

(43) 
(44) 

(45) 
(46) 



Thus we see that J may be interpreted as the contribution of this clump to the reciprocal 
magnification (relative to the 'empty' regions). Note that J may be negative; it vanishes on 
the caustics. Note also that 



J = 1 + L(A,A,) tr/C + L2(A,A,) det/C. 



(47) 



To find the average {A{Xs)), we have to mutiply (46) by the probability (28) of finding 
a clump at position (A, —x) and integrate over A and x: 



{A{X,)) = DliXs) 



dXno{l + zy / d^x{LtT}C + L^detlC) 



(48) 



Using (27), (33) and (36), we then find 



{A{X,)) = DliX,) 



3H^ng J dX{l + z)^L(A, A,) 



(49) 



which to this order is just (40). Note the crucial importance in deriving this result of the 
fact that the integral of the determinant of /C vanishes. 



5. The different averages 



To understand the important distinction between angular-size and area distances or 
magnifications, it is helpful to consider the shape of the surface of constant Zg (or constant 
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As, which as we argued is nearly the same thing) on the backward hght cone from the present. 
When only weak lensing occurs, it is obviously very close to being a sphere. When there are 
caustics, the light cone folds back and intersects itself. Then the world line of a comoving 
source will meet the cone several times, leading to multiple images. But it is important to 
realise that the different sheets of the surface are very close to each other; the time delay 
between the different images is always very small on a comsological scale. So the constant-^ig 
surface still lies very close to a sphere; it can essentially be described as a single sphere with 
some sections covered several times. We shall call this the source sphere. 

Now let us partition the celestial sphere around the observer into small pixels, labelled 
by n, with solid angles 5VLn (not necessarily all equal), small enough that the magnification 
may be taken to be constant within each. Of course Ylm = ^tt. For each n the hght 
arriving in the nth pixel originated from some small area 5Sn of the source sphere, with 

5Sn = An{\s)mn. (50) 

Note that some of these areas may be negative. 

Now, the result (40) established above shows that the total area of the source sphere is 
Y,5Sr, = Y,A^{K)5nr. = 47r(/l(A,)) = At^D^K). (51) 

n n 

With the demonstration that the area of the source sphere is the same as in the FRW 
universe of equal mean density, we provided a rigorous proof of Weinberg's argument, which 
was based on the principle of energy conservation (13). 

We can therefore define the solid angle Siln that would be subtended by SSn in the 
homogeneous FRW universe by 

Then the magnification relative to the FRW universe is 

" - MKY ^^^^ 

while the reciprocal magnification is 

The random- direction average of any quantity / defined on the celestial sphere is given 

by 

(/)rd= J]/n(A.)(5n„. (55) 
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In particular, of course, 



-) =^E^ = ^E^n = l. (56) 



Note that some of the terms in (56) are negative. If we consider instead the average of 
the modulus, we get 

\f^\ / rd n 

In fact, multiplied by Z)^(As) this gives the full area of the source sphere, counting the 
regions covered more than once with the appropriate multiplicity. We may define the ureal, 
as opposed to angular-size, distance 5area by 

DL.-D'{X,)(t^) ■ (58) 

\ l/^l / rd 

Obviously, -Darea(As) > D{Xs). 

We can also discuss random-source averages. To do that it is conveneient to choose the 
pixels in a special way. Let us first partition the source sphere into elements labelled by p, of 
area D'^{Xs)SQp. Here all SQp arc positive, and X^p^^p — ^tt. In general, each element may 
have several different images, though always an odd number, say 2rp + 1. Let us label these 
image pixels Sflpq, with q = 1, 2, ... , 2rp + 1. If we choose q so that positive and negative 
images alternate, then 5D,pi — 5D,p, SD,p2 = —SQp, . . . , and in general 

Mpg = {-y-^Mp. (59) 

It then follows that 

J]Mp, = 47r, (60) 

but 

The random-source average of some quantity / defined for each element of the source 
sphere is 

p 

Consider for example sources of given absolute magnitude randomly distributed on the 
source sphere, with UpSflp in element p. Then clearly 

Wrs = ^- (63) 
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What can one say about their average magnification? Actually, there are two possible 
meanings of the term. The first is appropriate if the various images of each source are 
unresolved. Then the total magnification of a source in Sflp is the sum of the magnifications 
of the separate images, each counted positively. 



It follows that ^ ^ 

(/^*°*)rs = — = — ^\Hpq\\5^lpq\. (65) 



47r ^ ^ 47r 

P,Q P,<1 



However, the signs of fj,pq and SQpq are the same, so we may remove the modulus signs, to 
obtain 

{l^'%s^^J2^^Pi-'^- (66) 

Thus from our earlier result, that the random- direction average of reciprocal magnification 
in unity, we have an important corolla ry, that the random- source average of total, direct 
magnification is also unity. Note, however, the important distinction, that every term in 
(65) is positive, whereas many terms in (56) are negative. 

The situation is different if the various images arc resolved (as may usually happen when 
the lenses arc clusters). Then the quantity of interest would be the average magnification 
of the individual images. Now clearly the total amount of radiant energy arriving at the 
observer is the same whether the images are resolved or not. The difference is that the total 
number of images is larger than the number of sources. To define the random-source average 
for resolved images, we should divide not by An but by the total solid angle the source sphere 
with all its folds subtend, namely AnDl^^jD'^ — 47r(|//|~^)rd- Hence 

Here again we have a relation between a random-source and a random-direction average. 

We conclude this section by estimating the magnitude of the difference between areal 
and angular-size distances for the single-clump approximation used in the preceding section. 
The reciprocal magnification If/j, when the light passes a single clump at (A, —x) is 

where J is given by (45). In going from the angular-size to the areal distance, according to 
(58), we have to replace /i by which in this case amounts to replacing J by | J|. Since 
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the integrand in (48) is just J — 1, this replacement yields 



^La=(l^(As)|)=i^^(A.) 



1 + j d\no{l + zY j (fx{\J\-l) 



(69) 



To be specific, we adopt essentially the modified isothermal sphere profile suggested by 
Brainerd et al. (2). In this distribution, the density at radial distance r is given by 

where a is the line-of-sight velocity dispersion, and s is a truncation scale, beyond which p 
falls rapidly, hke 1/r^. The total mass M of this distribution is finite, with 

M = (71) 

For our purposes it is more convenient to adopt a slightly modified version of the profile, 
designed to avoid the infinite density at r = 0. We take 

^^'^^ " 2ttG (r2 + a2)(r2 + s2)' ^^^^ 
where the core radius a <^ s. To be specific, we choose a = 0.01s. This profile shares most 
of the desirable properties of (70); over the wide region where a -C r <^ s, the density scales 
like l/r2, and the scattering angle ip can still be computed analytically. 

The results presented in the table below are for a flat universe with Hq = 70 km s~^ 
Mpc^^, Qm = 0.27 and fl\ = 0.73. We present illustrative parameters for galaxies with 
M = 1O"M0, a = 180 km s"^ and fig = 0.135, and for clusters with M = lO^^M©, a = 1000 
km s~^ and Qg — 0.0135, in both cases for sources at = 1 and Zg — 2. The table fists Dh 
and D, as well as the difference between ^area and D. 

The fact that the differences are much larger for galaxies than for clusters is due to 
the larger fraction of mass in the clumps, fig, and to the dependence on M and a. The 
relevant parameter is M/s^ or a^/M, which is larger for galaxies by an order of magnitude. 
At 2:5 = 1 the difference between areal and angular-size distances in not very significant, but 
by = 2, the effect of galaxies is to make the random-source average 3% greater than the 
random-direction one. 



6. Appropriate measures of magnification 

As we noted in the introduction, the appropriate measure of magnification depends on 
the type of source being observed, the mode of observation, and the manner in which the 
data are to be interpreted. 
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Table 1: The angular size and area distances {D and i^area) to a source at redshift Zg in 
a flm = 0.27, fl\ ~ 0.73 Universe where matter amounting to a normahzed density of VLg 
exists in isothermal spheres, each of total mass M and velocity dispersion cr, equations (71) 
through (73), and the rest of the matter is distributed homogeneously with a normalized 
density of fi^ — Q^. The quantity D represents the average reciprocal magnification over 
random directions in the sky, or the average magnification of randomly selected sources when 
the observed images for each source are unresolved. The parameter Dh^ which is the angular 
size distance under the scenario of the light path between a small emitter and the observer 
passing through only the homogeneous component (i.e. it misses all the clumps), is also 
tabulated. 



M 


a 




Zs 


Dh 


D 






(km/s) 






(Gpc) 


(Gpc) 


(Mpc) 




180 


0.135 


1 


1.74 


1.69 


13.1 


10" 


180 


0.135 


2 


1.96 


1.79 


62.6 


IQis 


1000 


0.0135 


1 


1.74 


1.69 


0.015 




1000 


0.0135 


2 


1.96 


1.79 


0.17 



Consider first a number count of discrete sources, say Type-lA supernovae. Clearly the 
expected number is proportional to the area of the source sphere. However, the number 
count can have two meanings. If the strong lenses are large — say galaxy clusters — so that 
the various images of a source are distinct and well resolved, then we may be interested in 
the total number of images, whether or not they come from the same source. In that case, 
wc need the area distance, Z)arca; the number of images is proportional to the total area of 
the source covered as we scan an area of sky. On the other hand, if the lenses are smaller 
and the various images are not resolved, we would be more interested in the total number 
of distinct sources. In that case, the appropriate measure is the angular-size distance. 

A different question wc may ask about discrete sources is their average magnification. 
Here we should use the random-source average. Again, there are two possible questions we 
might ask: What is the average magnification of each separate image? Or (if the images are 
not resolved) what is the average total magnification of the combined images? According 
to (66), the answer to the second question is exactly the same as in a homogeneous FRW 
universe of equaof equal mean density. On the other hand, the answer to the first question 
is obviously less, reduced by the factor (67). 

The situation with regard to a continuous source such as the cosmic microwave back- 
ground (CMB) is different. If an intervening lens magnifies a discrete source behind it, the 
source appears brighter, because it fills a larger area of the sky. This is not the case with 
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the CMB. The effect of magnification is that the radiation we see coming from one patch of 

sky originates from a smaller area of the source — in this case the last scattering surface. 
The surface brightness is unchanged. Lensing cannot change the temperature of a patch of 
sky, nor the temperature difference between two different spots. What it can change is the 
observed angular separation of these two spots. In other words, if there is hot spot on the 
last scattering surface, it may appear to us larger or smaller than it would in a homogeneous 
universe. 

For the CMB, the correlations between the temperatures at pairs of points arc controlled 
by the distance between source points rather than the area of source surface enclosed, so 
in considering for example the effect on the position of the acoustic peaks in the angular 
power spectrum we should use the angular-size distance. This is the measure for which 
the average is precisely the same as in a homogeneous universe of equal mean density. Our 
conclusion would not necessarily hold for very small scales, at or below the typical angular 
separation of lensed images, though such scales are below the limit of resolution of current 
CMB observations. 

This does not mean, however, that there could be no observable effect. The average 
magnification is unchanged, but the distribution around the average could be markedly 
affected with a separation between the mode and the mean. Lieu & Mittaz (9). Moreover, 
the acoustic peaks may be observably broadened, Lieu & Mittaz (10). One source of skcwness 
in the distribution operates on an angular scale such that the average region contains less 
than one clump. In that case, when we choose a random area of sky of the prescribed size, 
it is very likely to contain no clumps at all, in which case it will be demagnified. On the 
other hand in those rare cases when it does contain a clump, there will be a comparatively 
large magnification. So the peak may be shifted to a smaller scale. There has already been 
considerable discussion of this issue in the literature; see for example (Zaldarriaga & Seljak; 

1)- 



7. Discussion 

We have seen that so long as the weak-lensing condition is well satisfied, the average 
magnification or reciprocal magnification (in whatever sense) is essentially the same as in a 
homogeneous universe of equal mean (or pre-clumping) density, provided of course that our 
assumptions are correct, namely that the clumps are uncorrelated and well separated and 
that the gravitational effect of each can be treated in a Newtonian approximation. When 
the clumps are not well separated, then it may no longer be legitimate to assume that their 
contributions simply add, though we see no reason why that should not be true. We have 
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also ignored any time-dependence in the Newtonian potential, so the clumps must not be 
moving too rapidly, and must be small enough for the integrated Sachs-Wolfe effect to be 
negligible. Probably the most significant shortcoming of our approach is the assumption that 
the clumps are uncorrelated. Although we have considered explicitly only one size, there 
would be no real difficulty in including clumps of many sizes. But the effect of correlation 
between clumps of different sizes might well be important. 

When we approach the strong-lensing regime, on the other hand, the various averages 
are not the same, and in particular we must distinguish angular-size and areal distances. 
The random-direction angular-size average of the reciprocal magnification {1/ ii)^.^ is just the 
same as in a homogeneous universe of equal mean density, and as a corollary, the same is 
true of the random-source average of the total magnification of unresolved images, (/i*°*)rs- 
However, the random-source average of the magnification of resolved images, {11^^°^^^'^) rs may 
be significantly different, especially for the most distant sources. The difference between -Darea 
and D increases rapidly with increasing z. 



Appendix. Average of A 

In this appendix, we present a proof of the relation (40) between {A) and D. 

Let us consider a cylinder of fixed radius R centered on our fiducial geodesic, and 
extending out to A^. We choose R sufficiently large that clumps outside the cylinder have 
negligible effect; the value of R will drop out of the final answer. It is useful to consider in 
more detail the process of averaging over clumps. Clearly by (27), the probability that no 
clumps are found within the cylinder is 



Pq{\s) = exp 



Jo 



(73) 



The probability that exactly N clumps are found in the cylinder, at positions {{Xj, —Xj),j — 
1,2,...N} is 

N 

PoiKK Hi^ + Zj)^dXjd%. (74) 

Now what is the value of ^(A^) if clumps are present at these positions? Of course, 
if A^ = 0, then ^(A^) = -D/j(0, Xs)l- If there is just one clump, at (Ai, — Xi), its deviation 
effect is added, so we get 



A{X,) = ^^(0, A,)l + ^^(0, Ai)/C(fi)D^(Ai, A,). 



(75) 
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For N — 2, there are four terms, these two plus two more in which /C(x2) also appears. In 
general, we can write 

^(A,) = ^^,, (76) 
I 

where the sum is over all subsets / of {1, 2, . . . N}. Here, if / = {ii, ^2, . . . ip}, then 

p p 
MXs) = l[D„{X,^,X,.^,)l[jC{xi.), (77) 

j=0 j=l 

where Aj^ — and X^p_^_J^ — Xg. 

We can now recover our expression for (^(A^)). To do this we integrate (76) over the 
measure (74) and then sum over A^. Now when we integrate over Xi, then if i e /, by (33) 
we recover a factor of —AttGMI, whereas if i ^ /, we get simply ttR^. It is then convenient 
to sum separately over the number of elements P in the selected subset /, and the number 
in the complement, Q = N — P. It is easily seen that the sum over Q gives an exponential 
that precisely cancels the factor Po(As). Thus, again using (27), we obtain 

{A{Xs)) = iJ^i-oH'o^g) / dXj dX2... dXp 
p=o V J Jo Jxi Jxp-i 

p p 

n(i+^,)'n^'^(^^'Vi), (78) 

j=l j=0 

again with Aq = and Ap+i = Xg. It is easy to see that this is precisely the perturbation 
solution in powers of Qg of (38), so this again proves that (^(A^)) = -D(As)l. (Note that the 
required extra factor of (1 + Zj) comes from the initial condition (41) for Dh{Xj, Aj+i).) 

Next we turn to the computation of (A) = (det A) . Suppose as before that there arc 
clumps within our cylinder. Then from (76), using the identity det^ = |[(tr>l)^ — tr(>l^)] 
for 2 X 2 matrices, we have 

det^(A,) = ^ J] J][(tr^,)(tr^j) -tr(^,^j)], (79) 
/ J 

where / and J are arbitrary subsets of {1,2,..., N}, and Aj, Aj are given by (77). 

Now consider the effect of the integrations over Xj. As before, if i ^ / and i ^ J, the 
integral gives simply vri?^, and the sum over all such contributions generates an exponential 
that cancels Pq{Xs). Next, consider the indices that belong to one of the two subsets only, say 
i e / but i ^ J. Then )C{xi) appears only once in each term of (79), so we can immediately 
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perform the integral, and replace ]C{xi) with a factor — 47rGMl. Lastly, we are left with the 
indices, say {ii, 12, ■ ■ ■ is} that are common to both I and J. They contribute a factor 

^ J d^Xi^...(fxig ({tr[/C(fiJ . . ./C(fig)]}^ - tr{[/C(fjJ . . .A:(fis)]^}) 

= j d^Xi^ . . .d^Xigdei[K,{xi^) . . .'IC{xig)]. (80) 

But now in virtue of (36), this expression vanishes unless S — 

So finally we are left with a sum only over disjoint sets / and J. When wc have performed 
the integrations over a?,, all the matrices are proportional to 1, so the trace merely gives a 
factor of 2 which cancels the \ in (79). The remaining factors, for each of the two sums 
separately, are identical with those in (78), so we find 

(^(A,)) = [D{K)]\ (81) 

This concludes the proof. 
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